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ABSTRACT 

This study investigates the origin of interplanetary dust particles (IDPs) 
through the optical properties, albedo and spectral gradient, of zodiacal light. 

The optical properties were compared with those of potential parent bodies in 
the solar system, which include D-type (as analogue of cometary nuclei), C- 
type, S-type, X-type, and B-type asteroids. We applied Bayesian inference on 
the mixture model made from the distribution of these sources, and found that 
>90% of the interplanetary dust particles originate from comets (or its spec¬ 
tral analogues, D-type asteroids). Although some classes of asteroids (C-type 
and X-type) may make a moderate contribution, ordinary chondrite-like parti¬ 
cles from S-type asteroids occupy a negligible fraction of the interplanetary dust 
cloud complex. The overall optical properties of the zodiacal light were similar to 
those of chondritic porous IDPs, supporting the dominance of cometary particles 
in zodiacal cloud. 

Subject headings: zodiacal dust; interplanetary medium; comets: general; aster¬ 
oids: general; interstellar dust 


1. Introduction 


The purpose of this study is to investigate the origin of interplanetary dust particles 
(IDPs) taking into account the optical similarities and diversities between zodiacal light and 
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reflections from minor bodies in the solar system, snch as comets and asteroids. An enormons 
nnmber of IDPs are distribnted in interplanetary space. They are observable as scattered 
snnlight in the optical wavelength (zodiacal light) and as thermal radiation in the mid- and 
far-infrared wavelengths (zodiacal emission). The IDPs clond, occasionally referred to as 
a zodiacal clond, erodes on a time scale of 10^-10^ years (depending on the size and orbit, 
< 1/100 of the age of the solar system) dne to Poynting-Robertson drag, mntnal col l isions 
among the IDPs, a nd planetary pertnrbations flMa.nn et al.l l2006l: iGor’kavvi et al.l 119971: 


Dikarev et al. 


kg s ^ flGrnn et al 


20011). The mass-loss rate aronnd the Earth’s orbit is estimated to be ~ 10^ 


19851: iMann fc Gzechowskill2005l) . It is, therefore, natnral to think that 


ongoing dnst prodnction, snch as impacts or ice snblimation, is compensating for the erosion 
of the zodiacal clond. 

The origins of IDPs have been stndied throngh the spatial distribntion of the zodia¬ 
cal light. Early research for explaining t he spatial distribntion expected large contr i bntion 
from the asteroidal origin IDPs (see e.g. iDermott et al.l fjl99611 ). Later, iHahn et al.l (1200211 
compared the snrface brightness distribntion of zodiacal light, taken with the Glementine 
spacecraft onboard camera, to the inclination distribntions of comets and asteroids and sng- 
geste d that a significant fraction of dnst particles at 1 AU are of cometarv origin. Nesvorny 


et al. f 201ol ) farther performed a nnmerical simnlation for dnst particles ejected from six dif¬ 
ferent orbital gronps (asteroid families, main belt asteroids, Jnpiter family comets (JFGs), 
dormant JFGs, Halley-type comets, and Oort clond comets), and compared the brightness 
distribntion of the modeled zodiacal emission to that of zodiacal emission observed by an 
infrared space telescope. They snggested that 85%-95% of IDPs observable as zodiacal 
emission originate from JFGs. 

Althongh these recent stndies on the brightness distribntion favor cometary sonrces, 
little is known abou t the origin of IDPs in terms of their optical properties. Recently, 
Ishiguro et al.l (120131) derived the geometric albedo of IDPs by comparing the brightness of 
the Gegenschein, a par t of the zodiacal lig ht enhanced by backward scattering enhancement, 
to the infrared model (iKelsall et al.lll998l) . The research provided possibilities to study the 
origin of IDPs from a different aspect than previous studies. In this study, we considered 
the origin of IDPs through the comparison of the albedo and spectral gradient of zodiacal 
light with those of the potential parent bodies, and present a discussion based on previous 
studies. 
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2. Methodology 


The size distributio n of IDPs was studied through lunar microcrater counting and in 


situ flux measurements fiGrun et al.lIlQSSl : lDivinelll993l) . These studies suggested that the 


effective cross section of IDPs around the Earth’s orbit is dominated by large (10-100 /xm) 
particles. The opposition effect found in Gegenschein supports the idea that IDPs, which 
mak e UP the zodiacal light, are significantIv larger than the optical wavelength f ishiguro 


et al. 2OI3I : iBuffington et al.lboool) . Accordingly. we~could assume that the opticaF^operties 
of IDPs are similar to those of big objects, such as comets and asteroids. We thus postulated 
the albedo (A) and the spectral gradient {S') of the IDPs according to those of potential 
dust sources in the following discussion. 


2.1. Albedo and spectral gradient of zodiacal light 


The albedo of IDPs has been measured using several methods. iHapuerl fjl98011 com¬ 
pared th e zodiacal light br i ghtne ss to the IDPs model derived from the Lunar microcrater 
records. iLumme fc Bowelll (119851 1 derived the albedo value of IDPs, whic h can explain the 
polarization distribution of zodiacal light. iDumont fc Levasseur-RegourdI fll988l) compared 
the op tical and infrared brightness at the solar elongation of 90°. Recently, llshiguro et al. 


(120131) direct ly measured the ge ometric albedo of zodiacal light by comparing the optical 
and infrared (iKelsall et al.lll998l) flux at the anti-solar point. They deduced that the albedo 
of the smooth zodiacal light component is A=0.06 ± 0.01, after subtracting weak fine-scale 
features associated with asteroidal colli sional families, name ly, asteroidal dust bands. In this 
paper, we adopted the albedo value in llshiguro et al.l (1201311 . 


Small bodies in the solar system generally show linear spectra in a range around 4500- 
7500 A. It is useful to express the spectral index using the normalized reflectivity gradient, 

S' [% 10-3 A]: 

where S is the reflectance, defined as the flux density of an object divided by the flux 
density of the Sun at the wavelength A, and S and dS/dX denote the average reflectance and 
spectral gradient in the wavelength range, respectively. The spectra of zodiac al light have 


1996: 

Fixsen & Mather 

2002: 

Reach et al. 

2003: 

Ootsubo et al. 

2 OO 9 I: 

Tsumura et al. 


et al. 


2 OIOI ). Since there is no spectrographic data available in the optical wavelength range, we 
derived the optical spectral gradient S' by a log-line ar fitting using coni piled photometric 
data taken at different wayelengths around 4600 A ( Leinert et ahl 119981 1. The regression 
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formula for the ratio between the solar spectrum and the zodiacal light spectrum is given 

by, 

A 


h oc 


1.0 + 0.9 X log 


la 


0 


( 2 ) 


,5000 A 

where Ix and le denote the flux densities of zodiacal light and the Sun at wavelength A, 
respectively. Eq. (|2]) is applicable for zodi acal light in the spectral range of A < 5000 A at a 
solar elongation > 90° flLeinert et al.lll998h . From Eq. (|2]), we obtained the spectral gradient 
of IDEs as S'=8.5 ± 1.0 % • 1000 A ^ at 4600 A. Note that we derived S' at 4600 A in order 
to match the measured wavelength of the albedo fllshiguro et ah 2013 1. 


2.2. Data sources 


Turning now to the IDEs sources, we assumed that they origin ate from asteroid s and 
comets. In addition, some IDEs may originate from interstellar space flHahn et ahll2002l) . For 
asteroids, we considere d hve major taxonomic types, namely C-, X-, S-, B- and D-types 
flDeMeo &: Carrvll2013l) . as input data. Since the optical properties of cometary nuclei are 
similar to those of D-type asteroids (one taxonomic type of asteroids), we do not discriminate 
D-type asteroids from cometary nuclei. We thus assumed that IDEs consist of dust particles 
from six populations: C-type, X-type, S-type, and B-type asteroids as representatives of 
asteroids, cometary nuclei (including D-type asteroids), and interstellar dust. 

To create a template of the optical properties of six potential dust sources, we made 
use of catalogs of albedos and spectra of asteroids and comets. For aste roids, we used th e 
Asteroid Catalog Using AKARI(AcuA) catalog as a datase t of albedos flUsui et al.l 1201111. 


and the SMASSII catalog as datasets of spectral gradients f)Buslll999l: IBus &: Binzelll2002b 


Binzel et al.l 1200411 . We found 274 C-type, 222 S-type, 191 X-type, 40 B-type and 33 D-type 
asteroids archived in both catalogs. For C-type, B-type and D-type asteroids, which show 
no obvi ous absorption, we used the spectral gr adient values measured between 4350 A and 
9250A ( Bus fc Binzel 2002b; Binzel et ah 2004 1. For the S-type and X-type asteroids, w hich 


may have an absorption band around > 7000 A, we used the data at 4400 A — 7000 A (JB 


3 us 


19991 ). Albedos and spectral gradients of 10 cometary nuclei were compiled from various 
previous studies, shown in Table 1. For interstellar dust, we used the optical properties of 
average galac tic dust particles at 4600 A, that is, A=0.67 and S'=—23 ± 1 % • 1000 A 
flDrainel 120031) . We ignored some taxonomic types of asteroids, such as K-type, L-type and 
V-type as discussed in Section 4. 


Figure 1 shows the relationship between the albedos and spectral gradients for the 
IDEs and the potential parent bodies described above. In the diagram, the datum of IDEs is 














































5 


Table 1. Optical properties of cometary nuclei. 


Name 

type 

albedo 

spectral gradient 

references 

IP/Halley 

Halley type 

0.043 

7.5 

a, c, g 

2P/Encke 

Encke type 

0.050 

11 

h, j, r 

9P/Tempel 1 

Jupiter family 

0.056 

12.5 

q 

lOP/Tempel 2 

Jupiter family 

0.022 

20 

e, f 

28P/Neujmin 1 

Jupiter family 

0.025 

11.8 

b, k, p 

49 P/ Arend-Rigaux 

Jupiter family 

0.028 

10.4 

d, i 

67P/ Churyumov-Gerasimenko 

Jupiter family 

0.047 

10 

o, s, t, u 

103P/Hartley 2 

Jupiter family 

0.045 

8.1 

V 

162P/Siding Spring 

Jupiter family 

0.034 

9.2 

m, n 

C/2001 OGios (LONEOS) 

Halley type 

0.040 

9 

1 


Millis et al. (1988), e - 


Note. —a -ISagdeev et al.l fll986ll. 


A’Hearn et al. 


fll989l ) , h - ILuu fc JewittI fll990h 


fc Fer nan dez (1200211. 1 


J2nn6h 


Lamv et al.l (1200 


1 - 




LnnI 1119931 


b - Campins et al. ( 19871), c - Keller et al. ( 1987), d - 


fl989h. f - Jewitt fc Luu 1 1989k 


Abell et al.l fePObh. m -ICampins et al.l (12006), n - iFernandez et al. 


Campins et al.l (1200 


Li et al. 




Thomas &: Keller 


Fernandez et al.l (200011. k - Campins 




r - Boehnhardt 


et al. (2008), s - Lamv et al.l ( 20081) . t - ^biana et al.l ( 20081) . u - Kellev et al. ( 2009 ). v 
Li et al. 1 2013 ') 








































































































located within the population of comets and its spectral analog, D-type asteroids, suggesting 
that the major constituents are of cometary origin. In the following section, we will further 
investigate the contribution of each population using statistical analysis. 


2.3. Bayesian analysis 

When we chose objects from a type of population and calculated the correlation coeffi¬ 
cients between the albedos and spectral gradients, the absolute values were as low as -0.14, 
-0.04, 0.16, -0.36 and -0.11, for the C-type, S-type, X-type, B-type and D-type asteroids, 
respectively. Therefore, we considered these two properties, the albedo and spectral gra¬ 
dients, as being independent of each other. Within a population, we simply assumed that 
the albedo follows a log-normal distribution whereas the spectral gradient has a Gaussian 
distribution. We £t the distributions shown in Figure 1 to the model distributions shown in 
Figures 2 and 3. Because interstellar dust particles, which have limited contribution to IDPs, 
have optical properties that are very different from those of solar system objects, we assumed 
that the interstellar dust particles have a hxed albedo value oi A = 0.673 and a spectral 
gradient of S' = —23.2 ± 0.8 % • 1000 A , and did not consider its statistical distributions 
in the following analysis. 

We made a 2% grid for the possible combinations of fractional contributions from the 
source populations. Then, through linearly combining the probability distributions of the 
source populations according to the given contribution, we could generate a mixture prob¬ 
ability distribution for an optical property of a single dust particle. Different populations 
were weighted according to the average albedos. From the probability distribution for a sin¬ 
gle particle, we calculated the expected average values of both albedo and spectral gradient 
for the IDPs complex using Monte Carlo (MC) simulations. At the every grid point, a MC 
simulation with 500 sample particles was generated 3000 times. Under these conditions, 
the expected average values follow a Gaussian distribution, with a standard deviation of less 
than 10% of the uncertainty on the zodiacal light measurement. The probability of obtaining 
the measured optical properties was calculated at each grid point. By applying Bayesian 
inference with a flat prior, the probability was accepted as that of the assumed contribution 
of the grid point represent real situation. 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Albedo 


Fig. 1.— Spectral gradients S' with respect to albedos A of asteroids, comets and zodiacal 
light. Uncertainties of albedos are appended in the plot. The la measnrement nncertainties 
in spectral gradients are ordinarily abont 0.7 % ■ 1000 A . 

















































































0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Albedo 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Albedo 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Albedo 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Albedo 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Albedo 


Fig. 2.— Albedo di stribution of C-ty pe, S-type, X-type, B-type, cometary nuclei (including 
D-type asteroids) fjUsui et all 1201 ll) . Black solid lines are histograms for the given types. 
Red dashed lines are the lognormal distributions calculated from the mean and standard 
deviation of logarithms. 
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Spectral Gradient [% /1000 Angstrom] 



-30 -20 -10 0 10 20 30 

Spectral Gradient [%/1000Angstrom] 



-30 -20 -10 0 10 20 30 

Spectral Gradient [%/lOOOAngstrom] 


S-type Asteroids 




-30 -20 -10 0 10 20 30 

Spectral Gradient [% /I OOOAngstrom] 


Fig. 3.— Spectral gradient distrib ution of C-type, S-type, X-type , B-type, cometary nu¬ 
clei (including D-type asteroids) flBusI Il999l: iBns fc Binzell I2nn2bl ). Black solid lines are 
histograms for the given types. Red dashed lines are the Gaussian distributions calculated 
from the mean and standard deviation. 
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3. Results 

Table 2 shows the resulting contributions from the individual sources to the IDPs cloud. 
To derive the ranges (which are shown as plus and minus signs in the Table 2), we made 
contours with the same probability in the six dimensional grid, calculated the total proba¬ 
bility within contours around the most probable case, and derived the range with the 68.3 
% confidence interval. We found that cometary nuclei (including D-type asteroids) are the 
primary contributors (~ 94%) to the IDPs cloud as predicted in the Section 2.2. The re¬ 
maining part (~ 6%) is originated from the C-type, X-type and B-type asteroids. S-type 
asteroids and interstellar dust have an insignificant contribution to the IDPs (~ 0%). Figure 
4 shows the marginalized probability distributions of the four major populations. 


4. Discussion 


4.1. Feasibility of the method 


To assess the feasibility of our approaches described above, let us discuss the three 
following points. 

Firstly, we should consider the validity of the source populations. There are a wide 
variety of objects in the solar system, however, only six types of sources (five types of 
asteroids, comets and interstellar dust particles) are considered in this paper. Recently, the 
mass fraction from differe nt taxonomic typ e s of a steroids was studied using new multi-hlter 
photometric survey data. DeMeo fc Carry ( 20131) suggested that C-type asteroids account 
for more than 50% of the mass in the main belt. Although S-types, P-types, B-types, and 


Table 2. Contribution of 6 source populations to the IDPs 


Population 

contribution 

Cometary nuclei 


B-type asteroids 

4+l«% 

X-type asteroids 

21 ^ 4 % 

C-type asteroids 


S-type asteroids 


Interstellar dust 

01:1% 













Probability within 2% bin Probability within 2% bin 
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0.0 0.2 0.4 0.6 0.8 1.0 

Cometary Fraction 



0.0 0.2 0.4 0.6 0.8 1.0 

B-type Fraction 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 


X-type Fraction C-type Fraction 

Fig. 4.— (a) - (d) Marginalized probability for the fractions about cometary (D-type as¬ 
teroids), B-type, X-type and C-type asteroids. The probability for the vertical axis is the 
values integrated over 2% bin. 
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V-types have moderate fractions 10% each of the total mass of asteroids in the main belt), 
other asteroids snch as K-types, L-types and A-types only have minor contribntions of < 1%. 
P-types are inclnded in X-types in onr assnmption. Thns, we considered all bnt one, namely 
V-types, major asteroids in this paper. We conjectnre that V-types cannot co ntribnte to the 
IDPs clond becanse they have very large albedos (A = 0.30. lUsni et al.l (120111 )). In addition, 
the mass fraction of V-types is very small (0.01%) when we exclnde the largest objects in 
the taxonomic type (i.e. (4) Vesta). Meanwhile, the photopolarimeters on the Pioneer 10 
and 11 space crafts revealed that the zodiacal light brightness is negligible beyond 3.3 AU 
flTollerlllOSlI ). Those observations snggest that the contribntion from onter objects snch as 
the Kniper-belt objects (KBOs) may not be as large as those from asteroids, when we think 
abont the previons dynamic stndies point ing ont that the dnst particles from KBOs have 


peak densities ontside of the Jovian orbits flPoppe fc Horanvill2012l:IVitense et al.ll2014f) . The 
optical properties of the Centanrs show bimodality. Inactive Centanrs show nltrared spectra 
similar to KBOs, and active Centanrs have colo rs and albedos similar to cometary nnclei 
fjStansberrv et al.ll2008l: iMelita fc Licandroll2012|) . In this paper, inactive nltrared Centanrs 
were ignored along with KBOs and active Centanrs were treated as cometary nnclei. Some 
cometary nnclei have optical properties that are different from those of D-type asteroids, as 
in the cases of 95P/Chiron and 107P/Wilson-Harrington. We did not inclnde these kinds of 
objects in this stndy. We do not know how many snch objects exist, bnt these non D-type 
asteroidal nnclei are similar to other kinds of asteroidal gronps in terms of optical properties. 
Therefore each popnlation of an asteroidal gronp shonld be nnderstood to inclnde possible 
cometary nnclei whose optical properties are similar to the gronp. If we snbdivide the X- 
type asteroids into E-type, M-type and P-type, the resnlts remains same, with only the 
conhdence interval become worsen, becanse the optical properties of P-type asteroids are 
similar to those of D-type asteroids. 

Secondly, we shonld consider the time-evolntion of the optical properties via space 
weathering. We assnmed that the optical properties of dust particles resemble those of the 
source objects. However, it may not be true in some populations. Since the Poynting- 
Robertson lifetime of silicaceous dust p articles of 1 mm size is about 2 x 10^ years when 
released into circular orbit from 2.5AU fiMann et ahll2006l) while the time scale of the space 


weath ering is more than an order of magnitude shorter than the lifetime fjShestopalov et ah 
20131. ~ 7 X 10^ years for S-type asteroids.), it is reasonable to assume that surfaces on 


both silicaceous IDPs and S-type asteroids are altered by space weathering and therefore 
have similar optical properties. However, the space weathering of C- type, X-type and B - 
type asteroids are not well known, although there are studies, e.g. iMoroz et al.l (120041) . 
Therefore, we cannot clearly discuss the optical surface maturation of IDPs originated from 
these asteroids. Furthermore, cometary dust particles remain within the interplanetary space 
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longer than the active lifespan of cometary nuclei flLevison fc Duncanlll997l . ~12,000 years for 
the ecliptic comets), therefore the relation between the optical properties of cometary dust 
particles and the surfaces of active cometary nuclei is not direct. If we regard the cometary 
nuclei and D-type asteroids as identical, there are studies that imply that the spectra of D- 
type asteroids w ould not change sig nihcantly over time. D-type asteroids were found in the 
inner main belt (IDeMeo et al.l 120141) . and Phobos, possibly captured D-type asteroids, have 
optical p roperties of a D-type aste roid after remaining in the inner solar system for billions 
of years (jPaiola et al.ll2013l. l2014j) . Even though these objects have albedo values slightly 
higher than the average of D-type asteroids, their albedos and spectral gradients are still in 
the range of D-type asteroids. In other direction, according to the laboratory experiments on 
the Targish Lake meteorites, which have spectra similar to D-type asteroids, the continuum 
spectr um changed to the bluer direction after being exposed to laser radiation i Hiroi & 


Sasaki 120121 ). If these results can be applied in our case, the contribution of cometary nuclei 


would increase. If we think about the even now dominant cometary contribution, we can 
conclude that this assumption does not alter the conclusion of this paper. 

Thirdly, we should consider the effects of simplification in this study. We assumed that 
the optical proper ties are randomly dispersed within a population, however this may not be 
true. As shown in lUsui et al.l (120131) . there is a relation between orbital elements and optical 
properties. We want to emphasize that the differences in the optical properties between 
different types of sources are an order of magnitude bigger than the differences between 
sub-groups of different types of sources, as shown in Figure 1. Furthermore, we ignored the 
weak correlation between the albedo and spectral gradient. The relation was non-negligible 
for B-type asteroids. There is a possibility that the interstellar dust, which entered the solar 
system, has a different composition compared to the average dust particles in the Milky 
Way galaxy flMa.nnll2010l) . but we ignored this possibility. However, we want to justify these 
simplihcations because the contributions from B-type as teroids and interstellar dust are 
almost negligible. The SMASSII catalogue is not bias-free flMotliADiniz et al.ll2003l) . and we 
do not know the optical properties of unbiased populations, but we ignored the effect from 
bias. We hope that the large measurement uncertainties in the optical properties of zodiacal 
light cover the consequences of the bias. 


4.2. Comparison with IDP samples 


IDPs are nowadays collected in the Antarctic ice or the stratosphere around 20 -25 km al¬ 
titude using aircrafts, and they ar e well studied through laboratory investigations flBrownlee 


19851: lEngrand &: Maurettelll998l ). Because such particles should contribute to the zodiacal 
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light before they arrive on Earth, it is important to compare our result to IDEs samples. 
It is known that there are two major IDEs groups, which are referred to as ’’chondritic 
smooth” (CS) and ’’chondritic porous” (CE ). CS IDEs are composed of low porosity materi¬ 
als, predominantly hydrated layer silicates flSandford fc Walker! Il985l ). CE IDEs have large 
porosities of about ~70%. CE IDEs are dominated by anhydrous minerals. It is likely that 
CE IDEs origin ate from comets on the ground of mineralogical and petrographical properties 
(iBradlevI 120031 1. When the Earth passed through the dust stream of 26E/Grigg-Skjellerup 


(one of the JFCs), it was expected t hat 1-50% of th e total collected dust larger than 40 
/im could originate from that comet flMessengerl l2002l) . These were actually CE type IDEs 


with primitive ai ihydrous composition, supporting the assumption of that CE IDEs are of 
cometary origin flBusemann et al .1120091) . CS IDEs are thought to be derived from primitive 
(not differentiated) asteroi ds, because comets do not exhibit spectral signatures of hydrated 
silicate while asteroids do flMcAdam et al. 2015 b 


Bradley et al. (1199611 measured the reflectance spectra of IDEs samples in the optical 
wavelength and found that CS IDEs generally exhibit flat spectra with weak fall-off from 
6000 A towards 8000 A (similar to Cl and CM meteorites or C-type asteroids with S' ~ 0) 
whereas CE IDEs exhibit upward spectra {S' > 0) without remarkable curvature, although 
these IDE samples have a variety of albedo values and spectral slopes. Figure 5 compares the 
synthesized spectrum of zodiacal light based on our mixing mo del of small bodies w ith those 
of most typical CS IDEs (W7040A15) and CE IDEs (W030A5) fiBradlev et al. 1996 ). We also 
show the input spec trum in Figure 5, which was obta ined from a multi-ba nd photometry in 
Leinert et al.l (119981 ) and anti-solar point observation (jlshiguro et al.ll2013f) . The synthesized 
spectrum of zodiacal light is similar to that of CE IDEs (W030A5) in that it shows a low 
albedo value and positive slope {S' > 0) but it is different from that of CS IDEs in that 
it does not show a positive slope beyond ~ 6500A. It should be noted that the fall-offs 
below 4500 A in IDEs signals are artifacts o f the measurenients c aused by small size effects 
and should be ignored for the comparison (iBradlev et al.l 119961) . The spectral similarity 
leads to our assumption that the interplanetary dust complex is dominated by CE IDEs (i.e. 
dominance of cometary particles in zodiacal c loud). It is also curi ous to notice that CE IDEs 
tend to have cluster structures of 20-100 /im (IBradlev et al.lll996l) . The size is in accordance 
with th e effective size of zodiacal light dust particles evaluated by the IDEs size distribution 
model (jCrun et al.lll9851) . 


There seems to be a difference on the quantita tive estimates o f IDEs origins between 
laboratory investigations of IDEs samples and ours. iBradlevI (120031 ) studied 200 chondritic 
IDEs from the s tratosphere and found that about a half of them are classihed into CE 
IDEs. Similarly, iNoguchi et al.l (120151) investigated micrometeorite samples in Antarctica 
and suggested ~25% or even less are categorized into CE IDEs. Whereas we acknowledge 
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Fig. 5.— Comparisons of spectra between our synthesize d zodiacal light mod el (thick contin¬ 
uous line), CS and CP IDPs (dotted and dashed lines) (IBradlev et al.lll996l ). The observed 
reference spectrum of zodiacal light is also shown (see Section 2.1). Note the drop-offs in 
the IDP spectra at less than 4000 A artifacts. 
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that these laboratory investigations provide reliable results regarding the fraction of CS and 
CP IDPs fallen to Earth, we would draw attention to the sampling bias of the laboratory 
studies of IDPs. Asteroidal dust particles could be collected on Earth more selectively 
than cometary dust particles because of the orbital properties. The impact cross-section 
of asteroidal dust particles to the Earth is a few times larger than that of cometary ones. 
Furthermore, the impact cross-section of the Earth can be a few 1000 times l arger for dust 


parti cles trapped in quasi-satellite resonance, which favors asteroidal particles flKortenkamp 

boiah . 


4.3. Comparison with previous studies 


Our results are in agreement with kinematic, d ynamical studies based on spatial dis¬ 
tribution of zodiacal light. Numerical simulations of iNesvornv et al.l (1201011 concluded that 
> 90% of the zodiacal emission comes from JFC originated IDPs, and < 10% from the Oort 
cloud comets or asteroids. These results coincide with our results of > 90% cometary contri¬ 
bution. However, we could not optically distinguish different sub-populations of comets in 
this study. Therefore, we cannot know how large of a fraction of cometary IDPs originated 
from JFCs. 

We compared our results with infrared spectroscopic observations of zodiacal light. By 
using mean albedo values and contributio ns from this study and adopting a typical visual- 
infrared spectrum of the populations from iBus fc Binzell fl2002al) . we synthesized the model 
spectrum of zodiacal light, and extrapolated it to the near-infrared wavelength, shown in 
Figur e 6. At wavelengths shorter than 1.6 /xm, the observed NIR spectra of the zodiacal 
light (iTsumura et al.l 1201011 are similar to the synthesized spectra in large-scale, but have 
absorption-like dark wavelengths around 1.3-1.4 /xm. In this wavelength regime, the spectra 
of average B-type, C-type and X-type asteroids are bluer than those of zodiacal light, 
as expected in this study. The observed spectra have large uncertainties at wavelengths 
longer than 1.6 /xm. It is, therefore, hard to conclude clearly, but the observed spectra are 
closer to those of D-type, C-type or X-type asteroids than to those of S-type or B-type 
asteroids. Overall, the synthesized spectr um in this work is consistent with the obse rved 


spectrum from rocket-borne observati ons flMatsuura et al 


recent space observations o f AKARI (ITsumura et ahl 120131) . but it does not agree with the 


I 995 I: ITsumura et al.l l2010h and 


space observation of IRTS (jMatsumoto et al.lll996l ). The observed NIR spectra of zodiacal 
light are similar to those of D-type asteroids or cometary nuclei, which is in agreement with 
this work. Keeping in mind that we applied extrapolation through the model from optical 
observation, this match in the NIR reflected spectrum supports our conclusion of dominating 
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